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ABSTRACT 

Extrapolating from the solar system's asteroid belt, we propose that 
externally-contaminated white dwarfs without an infrared excess may be experi- 
encing continuous accretion of gas-phase material that ultimately is derived from 
the tidal destruction of multiple small asteroids. If this scenario is correct, then 
observations of metal-polluted white dwarfs may lead to determining the bulk 
elemental compositions of ensembles of extrasolar minor planets. 

Subject headings: planetary systems - white dwarfs 



1. INTRODUCTION 

An indirect but potentially powerful method to determine the elemental abundances of 
extrasolar minor bodies and planets is to measure the photospheric compositions of stars 
they pollute. Such studies in main-sequence stars are challenging because the host star's 
photosphere possesses intrinsic metals, making it difficult to determine the amount, if any, 
of contamination (see Gonzalez 2006). In contrast, heavy elements rapidly settle out of 
the atmospheres of white dwarfs cooler than 20,000 K (Koester & Wilken 2006), so that 
photospheric metals have an external source. We hope to identify which metal-possessing 
cool white dwarfs have accreted circumstellar rather than interstellar material in order to 
use these white dwarf atmospheres to measure the bulk elemental composition of extrasolar 
minor planets. 

There are strong arguments that the approximately ten single white dwarfs known to 
have an infrared excess (Kilic et al. 2006a, 2006b, Kilic & Redfield 2007, Mullally et al. 
2007, Jura et al. 2007b, von Hippel et al. 2007, Farihi et al. 2008a, b) are accreting tidally- 
disrupted asteroids. (1) The near-infrared excess discovered around G29-38 (Zuckerman & 
Becklin 1987) was shown to be produced by dust (Graham et al. 1990) that likely resulted 
from an orbitally-perturbed minor planet (Debes & Sigurdsson 2002) that passed within 
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the tidal-radius of the host star (Jura 2003). For other white dwarfs, the infrared excess is 
usually well modeled with a geomerically flat dust disk that lies within the tidal-radius of the 
star (Jura et al. 2007b). (2) The infrared spectrum of G29-38 displays strong 10 /zm silicate 
emission with a spectral shape that resembles dust emission from the zodiacal light rather 
than interstellar silicates (Reach et al. 2005). GD 362, the second white dwarf discovered 
to have a near-infrared excess (Becklin et al. 2005, Kilic et al. 2005), also displays a 10 /xm 
silicate feature characterisitic of asteroidal rather than interstellar dust (Jura et al. 2007a). 
(3) The white dwarfs with an infrared excess generally show when measured that F^(24 /im) 
<< F !/ (8 jum) (Reach et al. 2005, Jura et al. 2007b). There is no evidence for relatively cold 
dust grains outside the tidal-radius (Jura et al. 2007b). (4) The atmospheric composition of 
GD 362 is refractory-rich and volatile-poor (Zuckerman et al. 2007), reminiscent of material 
in the inner solar system and unlike the interstellar medium. While not as comprehensively 
studied as GD 362, there are at least three other white dwarfs which are makedly deficient in 
carbon relative to iron as is characteristic of asteroids and very unlike the Sun or interstellar 
matter (Jura 2006). (5) Asteroid belts similar the solar system's can supply the required 
pollution masses (Jura 2006). 

Only about 1/7 of all contaminated white dwarfs display an infrared excess (Kilic & 
Redfield 2007); the origin of the pollution of the more numerous contaminated white dwarfs 
without an infrared excess is uncertain. Although interstellar accretion could supply the 
metals in the atmospheres of contaminated white dwarfs (see, for example, Dupuis et al. 
1993, Koester & Wilken 2006, Zuckerman et al. 2003, Dufour et al. 2007), there are 
difficulties with this model. First, at least for helium-rich stars, the accretion of hydrogen 
relative to calcium must be suppressed by at least a factor of 100 (Dufour et al. 2007). 
None of the suggested mechanisms such as magnetic shielding that preferentially exclude 
hydrogen from the star have yet received any observational support (Friedrich, Jordan & 
Koester 2004). Furthermore, there is no correlation between metal contamination and the 
star's space motion as would be expected from Bondi-Hoyle accretion theory (Zuckerman 
et al. 2003, Koester & Wilken 2006). Finally, using a very simple model of grain inspiral 
because of Poynting-Robertson drag, Jura et al. (2007b) showed that the predicted fluxes 
at 24 jum for Bondi-Hoyle accretion can be more than 100 times greater than the observed 
upper limits. In view of these difficulties, the alternative of accretion of circumstellar matter 
merits consideration. 

Gaensicke et al. (2006, 2007) have identified two metal-contaminated single white dwarfs 
with gaseous dust disks that are likely the result of tidal-disruption of a large parent-body 
such as an asteroid. It is possible that in many cases, the dust orbiting a white dwarf is 
destroyed and therefore stars without an infrared excess may still experience accretion from 
tidally-disrupted minor bodies. In this paper, we elaborate upon this possibility. 



- 3- 



In contrast to hydrogen-rich white dwarfs where the gravitational settling time of metals 
is typically less than 100 years (see von Hippel & Thompson 2007), in helium-rich white 
dwarfs cooler than 20,000 K, the gravitational settling time may be considerably longer than 
10 5 yr (Paquette et al. 1986, Dupuis et al. 1993). Although very uncertain, below, we 
suggest that the lifetime of a dusty disk around a white dwarf might be ~10 5 yr. Therefore, 
at least for the helium-rich white dwarfs, the absence of an infrared excess might simply be 
a situation where a single tidally-disrupted asteroid has been fully accreted and the disk has 
dissipated but the contamination still lingers in the star's atmosphere. 

In §2 we present a qualitative overview of our proposed model where we suggest that 
contaminated white dwarfs without an infrared excess accrete from a gaseous disk created 
by the accumulation from many small asteroids. This contrasts with our model for the white 
dwarfs with an infrared excess which we expect to be the result of the destruction of a single 
large asteroid. In §3 we argue that asteroids larger than 1-10 km in radius are likely to 
survive the star's evolution through the Asymptotic Giant Branch (AGB) phase. In §4, we 
describe in somewhat more detail our model of white dwarf contamination by multiple tidal- 
disruptions of small asteroids. In §5 we discuss observational asssessment of this scenario. 
In §6 we discuss implications of this model and in §7 we present our conclusions. 



2. OVERVIEW 

Very little is known about extrasolar asteroid belts. Even the presence of warm dust at a 
few AU from the star does not tightly constrain the parent body populations (see Wyatt et al. 
2007). Here, assuming that many other main-sequence stars possess an asteroid belt similar 
to the solar system's, we present a single illustrative model to account for the white dwarf 
data. There are many unknown parameters; our calculations are aimed to establish there 
is at least one plausible scenario where multiple small asteroids can explain the presence of 
metals in white dwarf photospheres. If detailed abundance studies of polluted white dwarfs 
without an infrared excess are shown to be consistent with this hypothesis, then a more 
comprehensive investigation would be warranted. 

In our illustrative model, we suggest that the total mass of the asteroid belt when the 
star becomes a white dwarf is 10 25 g, somewhat larger than the solar system's asteroid belt 
current mass of 1.8 x 10 24 g (Binzel, Hanner & Steel 2000). We assume there is a steeply 
rising number of lower-mass asteroids, and we then compute that asteroids larger than ~1- 
10 km in radius plausibly survive the star's evolution through the asymptotic giant branch. 
Although we do not compute detailed models, we propose that the orbits of these asteroids 
are occasionally perturbed on a time scale of ~1 Gyr into orbits which lead them within 
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the tidal radius of the white dwarf (see Debes & Sigurdsson 2002, Duncan & Lissauer 1998) 
where the asteroid is rapidly shredded into smaller rocks and dust debris. 

Small asteroids are sufficiently common, that perhaps one is perturbed into an orbit 
where it is tidally-disrupted every 300 years. We argue that the accretion time on the white 
dwarf is long compared to this arrival time, so that white dwarfs with surviving asteroid 
belts inevitably possess a low-mass disk composed of disintegrated larger bodies. It is likely 
that newly perturbed asteroids arrive at a random nonzero inclination relative to the pre- 
existing disk. As the asteroid passes through the disk, solid material is rapidly sputtered 
into the gas-phase and a small asteroid is effectively vaporized. The matter gained by the 
disk from disrupted asteroids is lost by accretion into the atmosphere of the white dwarf, 
thus explaining the star's external pollution. In these systems, there is little dust and a 
negligible infrared excess. 

More rarely, a large asteroid is perturbed into an orbit where it is tidally disrupted. If 
sufficiently large, the destroyed asteroid dominates over the pre-existing gaseous disk. In 
these systems, a dust ring is formed and survives for perhaps 10 5 years; during this phase 
there is an infrared excess. 



3. ASTEROID SURVIVAL DURING THE STAR'S AGB PHASE 

Asteroids must survive the AGB phase if they are to pollute the star as a white dwarf. 
Rybicki & Denis (2001) discuss three processes that can lead to terrestrial planet destruction 
during the AGB phase. A planet can induce a tidal bulge in its host star which can then 
interact with the planet to drag it inwards. However, asteroids have such small masses, that 
this effect is not important. Below, we consider two other routes for asteroid destruction. 

3.1. Asteroid Destruction by Wind Drag 

Asteroids may spiral into the host star under the action of gas drag as they orbit through 
the wind of the AGB star. We assume that an asteroid is destroyed if it encounters its own 
mass in the wind. This destruction occurs because of hydrodynamic friction rather than 
sputtering which is negligible since the colliding gas atoms have kinetic energies of only a 
few eV (see Tielens et al. 1994). At distance D from the star in a spherically symmetric 
wind of mass loss rate M* with outflow speed, V win d, the gas density, p W i n d, is 



4nD 2 V wind 



(1) 
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As long as gravitational focusing is negligible, the rate of encountering mass, M enc , by the 
asteroid is: 

M enc = 71 R 2 ast PwindVast (2) 

where the asteroid of radius R ast has orbital speed V ast . For this evaluation, we assume that 
the asteroid's radius and mass are constant during the AGB phase. For circular motion: 

K„ - (^f (3) 

If the time scale for mass loss is slow compared to the asteroid's orbital period, then 
the specific angular momentum of the asteroid is approximately conserved and: 

DM* = DiMi (4) 

for an initial distance, D iy from a star of initial mass, Mj. As the star loses mass and M* 
diminishes, then the asteroid's orbit expands. Combining equations (l)-(4), then during the 
star's AGB phase, the total encountered mass, M enc is: 

M ,. -RL g 1/2 m? /2 



[■ R 2 ast G l/2 M-' 

, rnc — / M enc dt ~ ryr- (5) 

J lQV wmd Df /2 



In arriving at the last expression in equation (5), a simplified and slightly altered version of 
equation (3) of Duncan & Lissauer (1998), we assume Mj >> Mf where Mf is the final of 
the star (see, for example, Weidemann 2000). 

If p as t is the density of the asteroid whose mass is (4np ast Rl st )/3, then equation (5) 
leads to the condition that the minimum radius of a surviving asteroid, i? mm , is: 

R - 3M * /2gV2 (6) 

64fl pastV wind 

We only consider asteroids with initial distances from the host star greater than 2 AU 
because interior to this region, there may be a quasi-static region around the AGB star with 
a density much greater than given by equation (1) (Reid & Menten 1997). Beyond 2AU, 
we assume an outflow velocity of 2 km s _1 (Keady, Hall & Ridgway 1988), a speed that is 
appreciably smaller than the typical wind terminal speed of 15 km s _1 (Zuckerman 1980). 
We adopt p ast = 2.1 g cm" 3 , as inferred for Ceres (Michalak 2000). 

We show in Figure 1 the results for R m i n as a function of Di for initial stellar masses 
of 1.5 M , 3 M Q and 5 M Q . As can be seen from equation (6), the value of i? mm diminishes 
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rapidly as Di increases since less mass is encountered in the wind. Also, R m in increases with 
the intial mass of the star since the star ejects more mass for the asteroid to encounter. 
Depending upon the initial distance from the host star and its mass, we find that asteroids 
larger than 1-10 km in radius are likely to survive the drag induced by the AGB wind. 



3.2. Asteroid destruction by thermal sublimation 

Subjected to the star's high AGB luminosity, an asteroid might thermally sublimate, 
and we now estimate the minimum size asteroid that survives this process. If the star's high 
luminosity phase has duration, tAGB, then, the minimum size of a surviving asteroid is: 

Rmin = ^ tAGB (7) 

dR as t/dt is the rate at which an asteroid's radius shrinks because of sublimation: 

dRgst = o-(T) ^ 

dt Past 

where &(T) denotes the mass production rate from the dust per unit area (g cm" 2 s" 1 ) as a 
function of the asteroid temperature, T. Assuming pure olivine, then: 




a(T) = a \^e-W (9) 



where, converted to cgs units, & = 1.5 x 10 9 g cm" 2 s -1 and T = 65,300 K (Kimura et al. 
2002). Different materials may have different sublimation rates so our results are sensitive to 
the presumed asteroid composition. We adopt olivine as this appears to be the most likely 
carrier of the 10 fim silicate emission seen in G29-38 (Reach et al. 2005) and GD 362 (Jura 
et al. 2007a) as well as being common in the solar system. Assuming that the asteroids have 
a negligible albedo, we write that: 

T = ( ie^) " 4 < 10 > 

where L is the luminosity of the star and a SB is the Stephan-Boltzmann constant. While on 
the AGB, we assume the star has a luminosity near 10 4 L or larger for tAGB = 2 x 10 5 yr 
(Jura & Kleinmann 1992). 

Using equation (8)- (10), we show in Figure 2 the minimum size asteroid from equation 
(7) for stellar luminosities of 1 x 10 4 L and 2 x 10 4 L as a function of the asteroid's 
distance from the star. We neglect the orbital drift of the asteroid discussed above. Because 
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the thermal sublimation rate is temperature sensitive, asteroid survival decreases rapidly as 
the star's luminosity increases or for asteroids that orbit relatively close to the star. We see 
from Figure 2 that depending upon the luminosity of the star, asteroids with radii between 
1 and 10 km are likely to survive if they orbit beyond 3 or 4 AU. 



4. ACCRETION OF MANY SMALL ASTEROIDS 

Eventually, those asteroids that survive a star's AGB evolution, may have their orbits 
perturbed so that they come within the tidal radius of the white dwarf (Duncan & Lissauer 
1998, Debes & Sigurdsson 2002). When this occurs, the asteroid is shredded into dust. If 
the asteroid is relatively small and if there is a pre-existing disk from destruction of previous 
asteroids, then the shredded debris of the new asteroid is destroyed by sputtering in the 
pre-existing disk, and the circumstellar matter is largely gaseous. If, instead, the asteroid is 
more massive than the pre-existing disk, there is not enough gas to sputter the debris from 
the asteroid, and a disk with a large amount of dust at the inclination angle of the large 
asteroid is created. This more massive disk subsumes the low-mass pre-existing gaseous disk. 
Eventually, however, this large disk is accreted onto the star and the system relaxes to its 
steady state configuration of a largely gaseous disk. 

Denote the total asteroid belt total mass as M belt (t) : where t denotes the cooling time of 
the white dwarf. Let t or ut be the mean lifetime of an asteroid before it is perturbed into the 
tidal radius of the star. Assuming that all of the shredded asteroid ultimately is accreted, 
then the average metal-accretion rate, < M me tai >, is: 

<M metal >=^± (11) 

torbit 

If the circumstellar gaseous disk around a white dwarf has a lifetime, t gas , then the mass 
of a typical disk, M disk , is: 

M disk = < M metal > t gas (12) 

We propose that if an asteroid has a mass larger than M disk) then a single event dominates the 
circumstellar environment, and the star exhibits an infrared excess. If, however, an asteroid 
arrives with a mass lower than M disk , then this newly-shredded asteroid's dust debris is 
rapidly vaporized and the disk remains gaseous. 

We picture that all polluted white dwarfs possess in a steady state a gaseous disk, 
and, occasionally, a massive dusty disk. It is likely that a white dwarf takes relatively little 
time to enter its "normal" steady state with a gaseous disk. The first tidally-disrupted 
asteroid is shredded into dust which can persist for a long time since there is not likely to be 
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much viscosity in the dust particles. The second tidally-disrupted asteroid also is shredded 
into dust, but its orbital inclination is likely to be different from that of the first asteroid. 
Consequently, there will be grain-grain collisions approaching 1000 km s _1 within the tidal 
radius. These collisions are very destructive and large amounts of gas will be released. 
Eventually, as described below, the disk will evolve into a largely gaseous system. 



4.1. Duration of Gaseous Disks: t gas 

To progress further, we now estimate the lifetime of a white dwarf's circumstellar disk 
which we presume is controlled by viscous dissipation. We could extrapolate from cataclysmic 
variables to single white dwarfs. However, disks derive from destroyed asteroids probably 
have little hydrogen, and their compositions are very different from that of the disks in 
cataclysmic variables. For a gaseous circumstellar disk, the time for viscous dissipation is 
(King, Pringle & Livio 2007) 

tgas-D^ (13) 
av th 

where a is the usual dimensionless parameter, and V t h is ~/c_bT /m ast where ks is the Boltz- 
mann constant, T the temperature and m ast the mean atomic weight of the asteroidal matter. 
The most uncertain factor is a. As discussed by King et al. (2007), empirical estimates of a 
for thin, fully ionized disks are between 0.1 and 0.4, but theoretical simulations yield values 
of a that are at least an order of magnitude lower. In FU Ori stars, values of a as low as 
0.001 are inferred. The viscosity also may be a function of the gas ionization and decrease as 
the gas becomes more neutral. Here, because the disks around single white dwarfs may be 
long-lived, we assume that a is 0.001, but this value is extremely uncertain. We use models 
of the infrared emission from disks (Jura et al. 2007b) to adopt a characteristic disk size, 
D, of 3 x 10 10 cm. At this distance from the star, V as t = 5 x 10 7 cm s _1 , and at least for 
dusty disks, a representative temperature is 300 K (Jura 2003). If m ast equals that of silicon, 
then t gas ~ 5 x 10 4 yr. Although this estimate of the time scale is very uncertain, it is long 
compared to the time between arrivals of small asteroids within the tidal zone. 



4.2. Destruction timescale of Small Asteroids in a Gaseous Disk 

We now show that if a relatively small asteroid impacts a gaseous disk that the resulting 
infrared excess from dust debris has a short duration compared to the characteristic disk 
lifetime. In the solar system, asteroids in the main belt have a mean inclination, i, of 7.9° 
(Binzel et al. 2000), and in white dwarf systems, it is unlikely that a newly tidally-disrupted 
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asteroid has exactly the same orbital inclination as the pre-existing gaseous disk. As a result, 
as the debris from the asteroid orbits through the disk, it is eroded by sputtering. 

Let T, disk (g cm -2 ) denote the mean disk surface density: 

S*» = (14) 

If spherical grains of radius a move with speed v through a disk with number density, n, 
then the mass loss from sputtering is: 



— — = -y (ira nv) 15 

dt \ 3m ast J 



where y is the sputtering yield. If the grains and the disk have the same composition since 
they are both derived from asteroids, then each passsage through the disk erodes a layer of 
thickness, Aa pass , where: 

A(W = (16) 

4 p ast sin % 

With two passages per orbit, then in time t with t» t orbit , the minimum radius of a grain 
that survives, a m i n is: 

drain 2A(3p ass - (1^) 

torbit 

The infrared excess from a disk depends upon its optical depth which we now show 
rapidly becomes optically thin. Let x denote the opacity of the dust. Assuming that the 
grain optical cross section is given by its geometric cross section and making the conservative 
assumption to assume that all the dust particles have size a min , then: 

* = An ^ ^ 
^ Past (l m in 

If the mass of the asteroid is less than the mass of the disk, then the vertical optical depth 
through the disk, r 2 , is bounded so that: 

T z <X^disk (19) 

Using equations (16) - (19), then: 

3 sin i t oro it /or ,x 

Tz - — y — r (20) 

Light from the star of radius R*, illuminates the disk at a slant angle of R*/D, and the 
optical depth in the disk, r disk , to light at this angle is: 

D 3 D sin % t orbit 
r disk « r z — - — y — (21) 



-10- 



From equation (21), Tdisk diminishes with time as grains are destroyed. We now assume 
that the typical atom in the disk is silicon. The orbital speed at the disk is near 500 km s" 1 
and with an inclination angle of 7.9°, the characteristic speed and kinetic energy of a silicon 
atom impacting a grain are ~70 km s _1 and ~700 eV, respectively. Using the models of 
Tielens et al. (1994), we adopt y m 0.1. For D/R* = 30, as representative of a disk within 
the tidal zone of the star, equation (21) shows that r dist becomes less than unity within ~60 
orbital times. Since the orbital time of an asteroid might be ~10 years, this implies that 
there is appreciable infrared emission from the system for at most ~600 years. The true 
duration of an excess is probably much shorter since, for example, most incoming asteroids 
have masses appreciably smaller than that of the disk so that the bound given in equation 
(19) is conservative. Also, we have probably substantially overestimated the opacity of the 
dust debris by assuming all the particles have the minimum size. 



4.3. Duration of Dust Disks 

We propose that when a massive asteroid's orbit is perturbed so it is tidally destroyed, 
then a dust disk results. If analogous to Saturn's rings, this might only take ~100 orbital 
times or ~1000 years for the disk for form (Dones 1991). The subsequent evolution of the dust 
disk depends upon its viscosity which is highly uncertain. If the disk were composed purely 
of dust grains, then, like the rings of Saturn, the duration of the disk could be comparable to 
or longer than 10 8 yr. However, we anticipate that the dust disk is bombarded by orbitally- 
perturbed small asteroids and the resulting grain-grain collisions produce gas in the disk. 
For a mass distribution of asteroids given by equation (23) below, the median mass of a large 
asteroid that creates a dust disk is 2 M disk . A dust disk derived from the median-sized large 
asteroid accumulates by further collisions of small asteroids, its own mass in gas in a time 
scale of 2 t gas . Furthermore, the asteroid may have had a substantial amount of internal 
water which can be released into the gas phase as the asteroid is destroyed. Estimates for 
the percentage by mass of internal water in Ceres range between 0.17 and 0.27 (McCord & 
Sotin 2005, Thomas et al. 2005). If the viscous time scales inversely as the fraction of matter 
that is gaseous, then in the absence of adding any additional gas from small asteroids, tdust 
would be between a factor of 4 and 6 longer than t gas . Although extemely uncertain, we 
adopt on average that tdust ~ 3 t gas or perhaps 1.5 x 10 5 yr. 
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5. MODEL ASSESSMENT 

We now evaluate the model that asteroids are the main source of contamination many 
externally-polluted white dwarfs. In our "standard model" the mass of asteroid belt, Mb e i t (t) 
is taken to equal 10 25 g at the onset of the white dwarf phase when t = 0. We assume that 
there is a characteristic time, t orbit , to perturb an asteroid into a tidally-disrupted orbit. 
Therefore: 

M belt = MienfOje-'^' (22) 

This expression for an exponential decay of the asteroid belt is at best a rough approximation. 
Wyatt et al. (2007) discuss models where the mass varies as t" 1 because of mutual collisions 
while including orbital perturbations, more complex scaling models for the exponential decay 
may apply (Dobrovolskis, Alvarellos & Lissauer 2007). 

A critical parameter in our analysis is t or ut- To-date, simulations for the evolution of 
asteroid belts around white dwarfs have been somewhat limited, although Duncan & Lissauer 
(1998) and Debes & Sigurdsson (2002) computed various scenarios. For the future evolution 
of the solar system, Duncan & Lissauer (1998) followed the orbits of 5 asteroids, and found 
that 1 was always unstable. In two out of their four simulations, a second asteroid, Pallas, 
became unstable in 5 x 10 s yr. To reproduce the data, we adopt values of t or bu near 1 Gyr. 
However, a more simulations of for asteroid belt evolution during the white dwarf phase are 
needed. 

We assume that the asteroid belt has a mass distribution, n(M) dM, given by a power 
law so that: 

n{M)dM = A M~ 2 dM (23) 

where the constant A is derived by the normalization condition. This power law variation of 
M~ 2 is slightly steeper than the "classical" estimate for the solar system that n(M) varies 
as M" L83 (O'Brien & Greenberg 2005), but it is computationally convenient and adequate 
for systems about which we know very little. If M beU is the total mass of the asteroid belt 
with minimum mass, M min and maximum mass M max , then 

(24) 

As in the solar system, we assume the largest asteroid has an appreciable fraction of the 
total mass of the entire system, and we adopt M max (t) = 0.2 M be i t (t) Thus, initially, the 
largest asteroid has a mass of 2 x 10 24 g, somewhat larger than the mass of Ceres which is 
9.4 x 10 23 g (Michalak 2000). We assume that M m i n is independent of time and equals 2 x 
10 17 g, the mass of an asteroid with radius 3 km and density 2.1 g cm" 3 . In this model, the 
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total number of asteroids initially is 3 x 10 6 . With t or ut = 1 Gyr, this would result in an 
asteroid being tidally disrupted every 300 yr. 



5.1. Which contaminated white dwarfs have an infrared excess? 

In the context of our model, we now compute which white dwarfs may have an infrared 
excess. This requires estimating < M meta i > which we take as: 

< Mmetai > = 7^ e^'*** (25) 
^orbit 

To compare this quantity with observables, we estimate the white dwarf's cooling time from 
its effective temperature. We interpolate from the calculations of Winget et al. (1987) and 
Hansen (2004) and take 

U = 0.38 (|™) ' ,26) 

where tg denotes the cooling time (Gyr). This approach neglects the variation in stellar radii 
and masses. If a white dwarf exhibits an accretion rate above < M meta i >, then, according 
to our model, it is experiencing accretion from a single large asteroid and should possess an 
infrared excess. We show in Figure 3 a plot of accretion rate vs. stellar effective temperature 
for externally-polluted hydrogen-rich white dwarfs where the accretion rate can be reliably 
estimated. In this plot, we also show the value of < M meta i > for M = 10 25 g and t^t 
varying from 0.5 Gyr to 2.0 Gyr. The model and the data are consistent for most but not 
all stars. The adopted estimated mass in the asteroid belts of 10 25 g is larger than in the 
solar system; this point is discussed below in §6. 

While the data points in Fig. 3 mostly obey the predicted relationship, there are 
exceptions. The most deviant point is for WD 1633+433 with M meta i = 4 x 10 8 g s" 1 and 
T e ff = 6600 K (Koester & Wilken 2006). Perhaps this star has a particularly long-lived and 
massive asteroid belt. 



5.2. Relative numbers of contaminated white dwarfs with and without an 

infrared excess 

Any model to explain the contamination of white dwarfs should account for the relative 
numbers of stars with and without an infrared excess. We hypothesize in §4 that an infrared 
excess occurs when a tidally-disrupted asteroid is more massive than the ambient gaseous 
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disk. From equation (22), the number of such asteroids, N(M > M disk ), is: 

N(M > M disk ) « — — (27) 

The total rate at which these asteroids are perturbed into orbits where they are tidally- 
destroyed is N(M > M disk ) /t orbit . Therefore, if t dust is the duration of a dust disk, the 
fraction of time, fm, that the white dwarf exhibits an infrared excess is: 

f IR = N(M > M dtsk ) ^ (28) 

^orbit 

Using equations from above, this expression becomes: 

f IR = (tjt^A (in^-) 1 (29) 

Evaluation of expression (29) is insensitive to the exact values of M max and M min . From our 
estimates given above, we find that fiR = 0.062 (t dust /t disk ). In §4.3, we argued that t dust ~ 3 
t d isk- Although extremely uncertain, the prediction from equation (29) is therefore that Jm ~ 
0.19, consistent with the observational constraint determined by Kilic & Redfield (2007) that 
at least 14% of externally-polluted hydrogen-rich white-dwarfs possess an infrared excess. 



5.3. Composition of Accreted Material 

One motivation for studying the composition of contaminated white dwarfs is that it 
enables the measurement of the bulk elemental composition of extrasolar minor planets. For 
example, the photosphere of the white dwarf GD 362 was found to be markedly deficient 
in volatiles such as carbon and sodium relative to refractories such as calcium and iron 
(Zuckerman et al. 2007), a pattern reminiscent of the inner solar system where the Earth 
and asteroids also are very volatile deficient. 

Wolff et al. (2002) reported relative abundances of iron, magnesium, silicon, calcium 
and carbon in 10 helium-rich white dwarf. Although in many cases the uncertainties are very 
large, three of these stars exhibit a carbon abundance less than that of iron, and therefore 
a value of n(C)/n(Fe) that is at least a factor of 10 smaller than solar. These three stars 
have a composition similar to that of the Earth or chondrites; they do not appear to have 
accreted interstellar matter (Jura 2006). One of these three stars, GD 40, also displays 
an infrared excess (Jura et al. 2007b), consistent with the hypothesis that it is accreting 
chondritic material. However, the other two stars do not display an infrared excess (Farihi 
et al. 2008b) which, however, as discussed below, does not preclude the asteroid accretion 
scenario. 
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Desharnais et al. (2008) report FUSE data for 5 helium-rich white dwarfs. Of these 5 
stars, GD 378 has n(C)/n(Fe) = 2.5, within the errors not grossly different from the solar 
value of 8.3 (Lodders 2003). There are mostly only upper limits to the abundances for 
GD 233, a second star in this sample. However, G270-124, GD 61 and GD 408 all have 
values of n(C)/n(Fe) less than 0.3. These three stars are therefore candidates for having 
accreted chondritic material. G270-124 does not display an infrared excess shortward of 8 
/im (Mullally et al. 2007). As far as we know, GD 61 has not observed for an infrared excess. 
There is a hint of an excess at 7.9 /im for GD 408 (Mullally et al. 2007), but the signal to 
noise was not good enough for a definitive conclusion. 

Not all contaminated white dwarfs appear to have accreted chondritic material. G238- 
44 has been found to display n(C)/n(Fe) near 2.5 (Dupuis et al. 2007); this star does not 
display any infrared excess (Mullally et al. 2007). As with GD 378, the source of the external 
matter is unknown. 

There are at least two possible explanations for the absence of a near infrared excess 
around externally-polluted helium-rich white dwarfs. First, as discussed above, the cir- 
cumstellar disk could be largely gaseous. Second, the disk may have dissipated before the 
contaminants have settled out of the star's atmosphere. The settling time of heavy metals 
in the atmospheres of these stars with effective temperatures between 15,000 K and 21,000 
K ranges from ~3 x 10 4 yr to ~5 x 10 5 yr (Dupuis et al. 1993), comparable to or longer 
than our very rough estimate of the duration of a dust disk of 1.5 x 10 5 yr. 

5.4. Infrared excess from a largely gaseous disk? 

Above, we have adopted the view that a disk is either essentially all dust or all gas. 
However, it is possible that some dust survives in disks that are bombarded by numerous 
small asteroids. Consider the case where the interior of the disk is gaseous yet the outer 
portion remains dusty, and thus the disk distribution may have an especially large central 
hole. All but one of the ten white dwarfs currently known to have an infrared excess has dust 
as warm as 1000 K (see Jura et al. 2007b). This innermost temperature may be determined 
by the location where grain sublimation becomes rapid. However, G166-58 (also known as 
LHS 3007 or WD 1456+298), a metal polluted white dwarf with a metal accretion rate of 3 
x 10 7 g s _1 (Koester & Wilken 2006), was found by Farihi et al. (2008a) to have an excess 
at 5.7 /im and 7.9 /xm, but not at shorter wavelengths. Following Jura et al. (2007b), we 
compute a model with a flat disk to reproduce the data. We adopt a stellar temperature of 
7390 K and a ratio of the star's radius (9.1 x 10 8 cm) to the star's distance from the Sun 
(~30 pc), R*/D, of 1.0 x 10~ n . We fit the infrared data with an inner disk temperature of 
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400 K (corresponding to a physical distance from the star of 30 R*, Jura [2003]), an outer 
disk temperature of 300 K (corresponding to a physical distance from the star of 43 R*), 
and cosi = 0.4. This model disk lies well within the tidal zone which extends to ~ 100 R* 
(Davidsson 1999). The comparison between observations and the model is shown in Figure 
4; the fit is good enough that this model is a serious contender for explaining the data. Thus, 
this system may be in an intermediate phase of its disk evolution where the dust disk has 
an usually large hole because bombardment destroyed dust in the inner disk. 

5.5. Detectability of line emission from gaseous disks? 

Gaensicke et al. (2006, 2007) report emission lines from two disks; most polluted white 
dwarfs do not display emission lines even though we postulate the presence of a circumstellar 
reservoir from which they accrete. We suggest that the systems with detected emission are 
unusual because the gas has an unusually high excitation temperature, T ex . If the star has 
effective temperature, T*, and subtends solid angle, Q*, while the disk subtends solid angle, 
Qdisk, then a line is seen in emission if: 

B v (T ex ) fldisk ~ B V (T*) f2* (30) 

For SDSS 1228+1040, the analysis of Gaensicke et al. (2006) yields Vt disk /Vl* = 1.2 x 10 4 . 
Since T* = 22020 K, then for the calcium lines near 8500 A, the above criterion requires 
T ex > 1800 K. This temperature is much greater than expected for a passively heated disk 
material lying at 100 stellar radii (see Jura 2003), and some additional heating of the gas 
must be occuring. Gaensicke et al. (2006) show that the calcium lines are asymmetric and 
they explain this as the consequence of our viewing an eccentric disk. Because the orbits 
should circularize by mutual collisiions, it is plausible that the eccentricity is being driven 
by some external unseen planet. If so, this material also would be heated. Other disks may 
not be subjected to such extreme heating, and their circumstellar gas would be too cool to 
be detected. 



6. DISCUSSION 

There are a number of similarities between the asteroid belt in the solar system and the 
proposed asteroid belts around white dwarfs. We invoke asteroid belt masses of 10 25 g at 
the time the star enters its white dwarf phase, this mass is about a factor of 6 larger than 
the mass of the asteroid belt in the solar system. However, if stars form with a Salpeter 
mass function, n(M) dM, that varies as M~ 2 - 35 and if the lowest mass star that currently is 
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dying has a mass of 0.9 M , then the median main-sequence mass of a dying star is ~1.5 
M Q . There may be a linear scaling between the mass of a star's asteroid belt and its main- 
sequence mass. Furthermore, 1.5 M Q stars spend only ~3 Gyr on the main-sequence before 
becoming white dwarfs, significantly less than the current age of the solar system of 4.6 Gyr. 
Thus, when a 1.5 M star becomes a white dwarf, a simple linear extrapolation from the 
solar system suggests that it would possess 4 x 10 24 g, and thus approximately resemble the 
solar system. 

A minimum mass of an extrasolar asteroid can be estimated from the mass of contam- 
inants in the atmosphere of a contaminated helium-rich white dwarf. Since calculations for 
the mass in the outer mixing zone differ by as much as a factor of 10 (Dupuis et al. 1993, 
MacDonald et al. 1998), these estimates are uncertain. At the moment, the most polluted 
known He-rich white dwarf is HS 2253+8023 (Wolff et al. 2002), and the inferred mass of 
its accreted asteroid is ~10 24 g (Jura 2006), comparable to the mass of Ceres. 

The relative numbers of polluted white dwarfs with and without an infrared excess is 
explained by a mass distribution of extrasolar asteroids that varies as M~ 2 , approximately 
what is observed in the solar system. Thus extrasolar asteroid belts may possess a size 
distribution similar to that of asteroids in the solar system. 

Remarkably, in marked contrast to the Sun, carbon appears to be less abundant by 
number than iron or silicon in at least a few extrasolar asteroids. Similarly, as measured 
from chondrites (Lodders 2003), asteroids in the solar system have a similarly low carbon to 
iron abundance ratio. As argued previously (Jura 2006), determinations of the abundances 
in the atmospheres of white dwarfs can serve both to evaluate the model of interstellar 
accretion and, at least in some cases, to measure the composition of extrasolar asteroids. 

7. CONCLUSIONS 

We suggest that the contamination of many externally-polluted white dwarfs without 
an infrared excess results from tidal-disruptions of numerous small asteroids. This scenario 
is consistent with the hypothesis that the masses, size distributions and compositions of 
extrasolar asteroids are similar to those of the solar system's asteroids. 

This work has been partly supported by NASA. 
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Fig. 1. — Minimum surviving asteroid radius vs. initial distance from the star for stellar 
wind drag from equation (7). The solid, dotted and dashed lines are plotted for stars of 
initial masses of 1.5 M , 3.0 M and 5.0 M Q , respectively. 
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Fig. 2. — Minimum surviving asteroid radius vs. distance from the star for olivine sublima- 
tion (§3.2). The solid line and dashed line are for stars with luminosities of 1.0 x 10 4 L 
and 2.0 x 10 4 L Q , respectively. 
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Fig. 3. — Metal accretion rates vs. stellar effective temperature. The filled circles represent 
stars with an infrared excess while the open triangles represent stars without an infrared 
excess. The data are taken from Farihi et al. (2008a) and Jura et al. (2007b); the metal 
accretion rates are taken from Koester & Wilken (2006) and scaled according to the procedure 
in Jura et al. (2007b). The lines represent < M meta i > as a function of the star's effective 
temperature from §4.2. We assume initial asteroid belt masses of 10 25 g and the dotted, 
solid and dashed lines represent values of t dep i of 0.5 Gyr, 1.0 Gyr, and 2.0 Gyr, respectively. 
The solid circles are expected to lie above the curves while the open triangles are predicted 
to lie below the curves. 



-23- 



2.5 r — i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




Fig. 4. — Comparison of data showing la error bars (Farihi et al. 2008a) with the predicted 
stellar and disk emission for G166-58. The model fluxes from the star, the disk and the 
total are represented by the dotted, dashed and solid lines respectively. For clarity, the lower 
panel re-displays the dotted rectangle in the upper panel. 



